ftE  ET  PR^. 

/C^       y^   J^  NEW   YORK   UNIVERSITY 

<^   \jA  ]Y    ^  Institute  of  Mathematical  Sdences 

'^        vl/       X  Division  of  Electromagnetic  Research 


RESEARCH    REPORT   No.    BR-9 


A  Convergent  Asymptotic  Representation  for  Integrals 


J.  FRANKLIN  and  B.  FRIEDMAN 


CONTRACT    No.    AF  -  1  8(600)- 367 
DECEMBER    1  954 


VJ 


NEW  YORK  UNIVERSITY 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 


Research  Report  No.    BR-9 


A  CONVERGENT  ASYMPTOTIC  REPRESENTATION  FOR  INTEGRALS 


J.    Franklin  and  B.    Friedman 


IoaJ/    rkcLv^ k -^-^^  ct.B  ) 


Joel  Frariklin 


Bernard  Friedman 

-^/y      .  |5^/ ^E]^  YORK  UNIVERSITY 

-''^^-^       ^^"^-"-"^ )  pOURANT  INSTITUTE  -  LIBRARY 


Morris  Kline  _^,  . .        r     .      ^    ,  ,.,w 

„  X    T^-  X  251  Mercer  St.    New  York,  NY  '"'>"! 

Project  Director 


The  research  reported  in  this  document  has  been  made  possible 
through  support  and  sponsorship  extended  by  the  Office  of 
Scientific  Research,    H.  Q.  ,    Air  Research  and  Development  Com- 
mand,   U.S.A.F.,    Baltimore,    Maryland,    under  Contract  No.    AF- 
18(600)-367.     It  is  published  for  technical  information  only,    and 
does  not  necessarily  represent  recommendations  or  conclusions 
of  the  sponsoring  agency. 

New  York,    1954 


C.3. 


-  i  - 


Abstract 

This  report  presents  a  new  method  for  obtaining  an  asymp- 
totic representation  for  integrals  of  the  form  J    e"^  x^~  f (x)  dx 
when  p  is  large.  It  is  shown  that  if  f(x)  satisfies  certain  con- 
ditions this  representation  is  also  convergent.  Numerical  calcu- 
lations seem  to  show  that  the  first  term  of  the  representation 
gives  a  close  approximation  to  the  value  of  the  integral  for  a 
wide  range  of  values  of  p. 
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1,  Introduction 

In  many  applications  it  is  important  to  find  representations  for 
the  integral 

00 


/  e-P^  x°-^  f (x)  dx 


when  the  parameter  p  is  lsrge»  Here  c  is  a  complex  constant  such  that  Re  c  >  0  and 
f(x)  is  usually  an  analytic  function.  The  most  commonly  used  representation,  whose 
justification  is  due  to  Watson'-  -' ,  is  obtained  by  expanding  f(x)  in  a  power 
series  in  the  neighborhood  of  x  =  0  and  then  integrating  temi  by  term.  For  the 
sake  of  convenience,  we  shall  call  the  representation  so  obtained  Watson's  repre- 
sentation, even  though  it  has  been  known  from  much  earlier  times.  This  represent- 
ation has  several  disadvantages.  First,  in  many  cases  the  series  obtained  is 
divergent.  For  example,  if  we  take  f (x)  -  (1  +  x)"  ,  then  the  divergent  series 


I — I 

r  (c)p"°  Z.   (-)"  c(c+l)...(c+n-l)p- 


is  obtained.  Second,  the  representation  always  has  a  singularity  at  p  =  0  even 
though  the  function  defined  by  the  integral  may  be  a  regular  function  of  p,  as 
is  the  case  when  f(x)  -  (l  +  x)"  and  0  <  Re  c  <  1,  Third,  the  representation 
is  practically  useless  for  computation  when  p  is  not  very  large  but  only  in  the 
range  from  2  to  5* 

In  this  paper  we  present  a  method  for  expanding  the  above  integral  into 
a  series  which,  in  many  cases,  is  both  asymptotic  in  p  and  convergent  for  p  >  Oj 
consequently,  the  series  may  be  used  for  computation  for  all  non-negative  values 
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of  p.     Moreover,  for  moderate  values  of  p  the  first  term  alone  usually  gives 
more  accurate  results  than  the  asymptotic  representation  obtained  by  Watson's 
method.     A  further  advantage  of  the  first  term  is  that  it  is  not  singular  for 
p  -  0  if  f(x)    goes    to  zero  like  x~    when  x  goes  to  infinity.    However,  it  should 
be  remarked  that  the  series  obtained  in  this  paper  is  much  more  complicated  than 
that  obtained  by  Watson's  method  and  that  in  many  cases  the  complete  series  can 
be  obtained  only  with  great  difficulty. 

2*     The  First  Term  Approximation 

Let  us  denote  the  integral 

-00 


/ 


^-PX      ^C-1   ^(^)^ 

o 

by  the  symbol  L   [f (x)] .     If  we  wish  to  make  explicit  the  fact  that  this  integral 
c 

is  a  function  of  p,  we  shall  use  the  symbol  L   [f(x)j  pj» 

As  is  well  known,   the  first  teim  of  Watson's  representation  for  L   [f(x)] 
is  f(0)  n(c)p~'^.     The  error   incurred    in  using  this  tenn  for  the  integral  L  (f(x)] 

^^00  00  1 

E  -  /     e-P^  x*^^[f(x)-f(0)]dx  -  /     e-P^x^  dx     /     f '  (xv)  dv 

if  we  assume  f(x)  is  continuously  differentiable.     Suppose  that 

lf'(x)|   <M 
for  0  <  X  <  oo;  then 

|E|  <Mr  (c  +  1)  p-°"^      . 

This  shows  that  the  ratio  of  the  error  to  the  first  term  in  Watson's  representation 
is  of  the  order  p"  . 
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A  better  approximation  to  the  value  of  the  integral  L^[f(x)]  will  be 
given  by  the  expression 


(1) 

Now,  the  error  is 


f(|)  r  (c)  p"". 


00 


\  ' 


-px   c-1  J 


e  "^   X 


f  (x)  -  f  F)  , 

K.  P 


An  integration  by  psrts  shows  that 


(2) 


/oo 


-px       c    _d_ 


e  ^       X 


dx 


f(x)  -  f^) 

L         ^-P  .: 


'1    /"^  e-P^  /  dx     /  vf"  [^  *  v(x  -  ^)]  dv. 


where  we  have  asstmed  that  f(x)  has  a  continuous  second  derivative.  If  we  assume 
also  that 

If"  U)l  <M^ 

for  0  <  X  <  00,  then  we  get 


O) 


ie^i  <M^r  (c  + 1)  p" 


-0-2 


This  estimate  shows  that  now  the  ratio  of  the  error  to  the  first  teim 
is  of  the  order  p"  ;  conseoiiftntly,  we  expect  that  (l)  will  give  a  closer  approxi- 
mation to  the  value  of  L^  [f (x)]  than  the  first  term  in  Watson's  representation. 
The  following  illustrations  show  the  accuracy  of  our  approximation 

1)  Suppose  that  f(x)  =  (1  +  x)   and  c  =  1}  then  from  (l)  we  get 

/oo 
e"^   (1  +  x)-^  dx  '^^  (1  +  p)"^. 


-  li- 


In  the  following  table  we  compare  for  real  and  imaginary  values  of  p 
the  value  of  the  integral,  which  we  denote  by  L^,  with  the  value  of  the  approxim- 
ation: 


p 

1 

2 

3 

h 

5 

h 

.5963 

.3613 

.2621 

.2063 

,noh 

app. 

.5000 

.3333 

.2500 

.2000 

.1667 

p 

i 

2i 

3i 

Ul 

5i 

h 

.3U3-.622i 

.212-.399i 

.079-.292i 

.050-.229i 

,03U-.l88i 

app. 

.500-.5ooi 

.200-.li00i 

,100-.300i 

.059-.235i 

.039".192i 

2)  Consider  the  case  where  f(x)  =  (l+x)  ''  and  c  =  1,  We  get 


00 


.-1/2 


(l+x)  ^  dx  ^  (p  +p) 


2..^-l/2 


In  the  following  table  we  compare  the  value  of  this  integral,  which  we 
denote  by  I„,  with  the  corresponding  value  of  the  approximation  *. 


p 

h 

5 

Ui 

51 

h 

.2263 

.18U2 

,0268-.2laiii 

.0179-. 19521 

app. 

.2236 

.1826 

.0301-.  2lMi 

.0195-.1971i 

3)  As  a  final  illustration,  we  shall  obtain  an  approximation  for  the* 
I2] 


Hankel  function.  We  have' 


,^i)(y) .  ^  y"  ,^y  -*?  d^  . 
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Put  cosh  V     =  1  +  xj  then 


«("(„  .  ^  L 


where 


00 


X--/2  e^y^  (2.x)-^/2  dx. 


-1/2 
If  we  put  c  -  1/2,  p  -  -iy,  f (x)  -  (2  +  x)  '  then  the  integral  L  is  exactly 

the  type  we  have  been  considering^  consequently,  from  (l), 


L-r(i)(.iyr^/^2.^r^/^ 


^tS^^ 


1/2 


1/2 

-  (2^    e 


i,A  y-l/2(^  ^  ^-l/^  ^^^   arctanCUy)-!] 


We  see  then  that 


H^^^(.) 


where 

F     -     fi  ♦  ^:-.)-l/^    exp[-i  arctan(Uyr-M    . 

The  usual  asymptotic  expression  for  H       (y)   is  obtained  by  taking  the  limit  of 
F  as  y  approaches   00,   that    is,  F  «  1« 

The  followinf^  table  shows  the  closeness  of  the  approxmation  when 
F  is  not  replaced  by  1. 


y 

0.6 

2 

u 

5 

i^"(y) 

.912-.309i 

.22U+.5lOi 

-.397-.017i 

-.178-.309i 

app. 

.9l8-.370i 

.228+,5Ui 

-.398-.017i 

-a78-.309i 
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3.  The  Complete  Expansion  of  L  [f  (x)] 

It  follows  from  Equation  (2),  which  gives  the  error  in  using  one  term 
for  L  £f(x)J,  that  this  integral  can  be  written  as 


(U)  Ljf(x)]  =  f(^)  r(c)  p""  *   p-^  I^c+l&i^^)  ^ 

where 

(5) 


-  r  f (x)-f (^) 

^l(-)  -  h   i  ^  r    . 


X  -  — 

p 


It  is  clear  that  in  (U)  we  may  replace  c  by  c  +  1  and  f(x)  by  f ,  (x)  to  get 


^^^      Vl^V^^]  "  fl(~i)  r  (c+l)  P"°"^+  p"^  L^2  ^2^^^] 


where 


fgCx) 


^  \  c+1 


x  - 


Combining  (U)  and  (6)  and  continuing  this  process,  we  find  that 


(7) 


L  [f(x)]  =  p-=  r  (c)  f" 


:(c+l)...(c+k-l)  p"*^  fj,  (—--) 


*  ""Vn  &„(=='] 


where 


(6)      f^(x)  .  f(x),  f^,^(x)  =  |j 


f'kt^Ktfi 


L 


X  - 


c+k 


If  we  let  n  go  to  infinity  in  (7),  we  obtain  the  infinite  series 


00 


(9) 


p"°  r(c)  r  c{c+i)...(c4k,i)  p'^\(~) . 

k»0 
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This  series  is  asymptotic  in  p  in  a  sense  which  is  rasde  precise  by  the  follo;d.ng 
theorem: 

Theorem  I: 

Put  c  -  a  ♦  i  p,  vrtiere  a   is  real  and  positive  and  p  is 

real.  For  P  j^  0,  suppose  that  f  (x)  is  analytic  for  Re(x)  >  0  and 

I       (2N) 
that  f,  f  !•••,  f     satisfy  the  inequalities 

(10)  |f^^\x)|  <  M  e^"'""',  hO,l,...,2N 

where  M  and  |i.  are  non-negative  constantst,  Then  for  n  =  1,2,.«.,  N, 
we  have  as  p  -*  oo 

N-1 

(11)  L   [f(x)}  =  p-'Rc)  L  c(c+l)...(o.k-l)p"^\(^)  +  0(p-2^^. 

k-0 

If  p  =  0,   that  is,   if  c  is  real   and  positive,    the  assmmpticn  that  f (x) 
is  analytic  for  Fte   (x)  >  0  may   be    replaced  by  the  assumption  that  f(x) 
is  2N  times  continuously  differentiable  for  x  >  0,  and  then  (11)  holds 
with  a  replaced  by  c. 

In  many  cases  the  series  (9)   is  not  only  asymptotic  but  also  convergent, 

We  shall  prove 

Theorgn  II; 

Suppose  f (z)  ■  f (x+iy)  can  be  represented  in  the  form 

f(z)  -  /   e"^"^  dT(t),  x>0, 

o 

where  Y(t)   is  a  complex-valued  function  which  is  of  bounded  variation 
in  each  finite  interval  0  <  t  <  T  and  which  satisfies  the  inequality 

|T(t)l   <M,  t  >  0 

Then  for  p  >  0  and  Re  c  >  0,  the  series  (9)  converges  to  L  [f (x)J , 
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This  theorem,  when  applied  to  the  integrals  used  as  illustrations 
in  the  previous  sections,  shows  that  the  series  converge. 


Proof  of  Theorem  I.  Equation  (?)  shows  that  we  need  prove  only 
Equation  (?)  shows  that  we  need  prove  only 


(12) 


^e+N  Tyx)]  -  0(p-«-^). 


First,  we  prove  by  induction  that 


(13) 


|f^^)(x)|  <Mj^e^'l^l 


for  h  =  0,  1,  2,,,,,  2N-2k  ^Jid  k  -  0,  1,  2,.o.,N.  Here 


M.  »  M  exp 


p-S 


k^*k(|c|-l) 


From  (8)  we  have 

1 

^k+l^""^  "1^/4  (p*"'"(c+k)(l-v)  +  vx]  dv 

(lii)  1 


y^f"  [p"-^(c+k)(l-v)  +  vx]  dv  . 


Suppose  now  that  (13)  holds  for  some  value  of  k;  then  differentiating  (14) 
k  times  with  respect  to  x  and  using  (13)  we  find  that 


f^^^  (x) 
1^1  ^^^ 


\  f    v*"*"^  exp  [n  I X I  +p"^(  I  c  l*k)]  dv 


|i|x| 


^\*1  ^ 


for  h  =  0,  1,  2,.,.,  2N-2k-2,     By  hypothesis,  this  inequality  holds  for  k  -  0, 
consequently  it  holds  for  0  <k   <  N, 
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Now,   using   (13),   we  see  that 


00 


lc.i.[V->]|-/«'''^'"^''N<''^ 


00 


o 

-  0  (p-'^). 


-px     N+a-1  |jx 


e  *     X  M^  e'^  dx 


irtiich  proves  the  theorem* 


Proof  of  Theorem  II .  Be'^ore  proving  this  theorem,  we  shall  obtain  a  more  conven- 
ient  representation  for  f.  ( )»  Suppose  that  fv.(x)  ■  e   j  then  from  (lU)  we 


have 


Vl^-^ 


or,  if  we  put  W  «  av. 


-  /  va  expl av(x  -  — ) 


dv 


(15) 


fl^+l(x)  -  /  W  exp 


-(c+k)(u-W)p' 


e    dW« 


-tx 
We  can  now  show  that  if  f  ■  e   ,  then 

o      * 


W 


(16) 


h-1 


f^(x)  -  //  •••/   W^W^" 


-W  X 

•W  A  e  "  dW  dW  ,**MW, 
n  n       n  n-1    1 


o  o 


for  n  ■  1,  2,...,  where 


A^-  exp  -ctp"-^  -cp"^ZI  \  *  (c+n-l)p"-Hj^ 


The  proof  is  by  induction.  Formula  (l5)  shows  that  (16)  holds  for  n  ■  1.  Suppose 


that  (16)  holds  for  n-1,  2, •••!€;  then,  using  (l$)  with  a  -  W  and  interchanging 


-  10  - 


the  order  of  integration  (this  is  justified  if  x  >  0),  we  get 

*     ^  Vl 


^k.l(^) 


W  •••>/     A  dW  •••dWj^ 


o     o 


\     r 

•  y    W  exp  - 


(c+k)(Wj^-W)p' 


-1 


e       dW, 


Clearly, 


k^  exp[-(c+k)(Wj^  -  W)  p-lj   -  Aj^^^  . 


This  shows  that  (16)  is  valid  for  n  »  k  +  1  and  thus  for  all  values  of  n. 
By  the  hypothesis  of  Theorem  II, 

00 

f(x)  -  /  e'*'''  dY(t). 

o 
To  obtain  f„(x).  '"^  note  first  that  we  may  write 

fn(x)  -  Pnfr(x)] 

where  P  is  a  linear  operator  acting  on  the  space  of  2n-times  diff erentiable  func- 
tions. In  fact,  the  operator  P  may  be  defined  for  n  »  1,2,,..,  recursively  from 


(111)  as  follows:             , 

r 

0 

f' 

and 

p"  (c+n-l)(l+v)+vx( 


dv 


P^[f(x)]  -  f(x). 

If  X  >  0,   it  is  easy  to  justify  the  fact  that 

-     00    _.  -V  00  _.  - 

^n^""^  "  ^n    y^  «"     d^(t)     -/dl(t)  Pn    e 

Lo  -^       o  L         . 


Using  formula  (16),   we  get 

00  t  W, 

dl(t) 

-1. 


Vl  -W  X 

•••  /        W  •••¥  A  e     "  dW  ♦••  dW-. 
/  Inn  n  1 


Put  X  "   (c  +  n)p     f  then  we  get 


(£> 


(18) 


f„(^)  -  /  dT(t)  e-^^/^  l(t), 
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where 

+     w  \ 

,n-l 


t     VL  W 


I(t)=/     /      •••/       W  •••  W^  exp(-  i  y[  Wj^)  dW  •••dW  , 


(19) 


The  region  of  integration  in  (19)  is  an  n-dimensional  pyramid  with  vertex  at  the 
oi^igin.  Since  the  integrand  is  a  symmetric  function  of  its  n  argiiments,  the  value 
of  the  integral  over  the  n-dimensional  cube  0  <  W  <  t,  (k  »  l,2...,n)  will  be 
exactly  nl times  the  value  of  I(t).  This  shows  that 


I(t)  .lr\    /we-^/PdW 


■)■ 


1   2n  J,   f.   ^  ts   -t/pl 
"ITT  P   jl-  (l--)e  ^P>   . 

Substituting  this  value  for  I(t)   in  (18),   we  obtain  the  desired  results 


(20) 


f„(^)  -  ^  P^y  *df  (t)  e-*/P  |l  -  (1  .  p  e-^/pj"  . 

With  the  help  of  this  representation  the  series  (9)  may  be  written  as 

^o 

where 

(22)  J?(t)  «  1  -  (1  +  t)  e'"*^  . 

To  prove  Theorem  II,  we  must  now  prove  that  this  series  is  convergent  and  that 
its  sum  has  the  value 

00  00  00 

(23)  Ljf(x)]  -  /  x^-VP^f (x)dx  -  /  x^-VP^dx/  e-^dT(t) 


'O 

00 


dl(t)  (1  +  i)"° 

p 
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The  function  0(t)  increases  steadily  from  zero  to  unity  as  t  increases 
from  0  to  CD,     From  the  binonial  theorem  it  follows  that 


(2U)      (1  *  i)-=  -  e-°*/P  h   -  0{p[    .  e-^*/P  ^  r  *  k  "  7  ^^^P 
for  0  <  t  <  00.  Here  we  have  used  the  fact  that 


Inserting  (2U)  into  (23)  and  integrating  the  first  n  terms,  we  find  for  the  in- 
tegral in  (23)  the  representation 


(25) 


-  00 

n«l 


r  C*r')/«-'*^'^«''(FW*«„ 


where 

(26)  R^  y  e-^^/P  r  ('  '  k  "  ^)  ?^''(|)dl(t). 

o 

From  (21)  and  (25)  it  follows  that  the  theorem  will  be  proved  when  it  has  been 

shown  that  R  -^  0  as  n  -!>  oo« 
n 

For  the  sake  of  convenience,  put 

co(t)  -  I(tp)j 

then  co(t)  is  of  the  bounded  variation  in  each  finite  t-interval,   and   |oo(t)|  <  M. 
Replacing  t  by  tp  in  (26),   we  get 

°°  00      /  \ 

(27)  «n-/     "*    H     (  '  *  ^  -  ')  /(t)  dco(t). 
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Wa  will  now  prepare  to  integrate  by  parts.  The  series  in  the  integrand  may  be 
differentiated  term  by  term;  the  derivative  of  the  integrand  is,  therefore,  equal 
to 

(28)     H- E  ("  T  ')  ^'(t)  ^  E  k  (°  *  ^  -  ')  <^\t)   te-*l  . 


Using  the  identity 


te 


-t 


T^ 


l-0(t)>. 


we  may  write  the  second  summation  in  (28)  in  the  form 


(29) 


i^|<=*r^)[- 


(t)  -  fit)) . 


After  partial  summation,  (29)  becomes 


t 

1+t 


n 


;*n-^ 


k»n 


(- ^'(^ :  i) -4' *  r  l^"'" 


(30) 


t  .  nh  * ;  -  M  rht)  *  gj 


00 


1+t 


■n 


"  *  IJ  - 1  At) . 


Here  we  have  used  the  identities 


y{i:'i]--(''i-']-  f  *r^) 


(k+  1] 
Using  (30),   we  may  now  write  the  derivative  (28)  in  the  fonn 


(31) 


te 


-ct 


1+t 


•  n  I  °  *  '^  -  1)  0^1(t)  .S^    f^[<^*l'^]  0\t), 


n 


1+t 


k-n 


-  lU  - 

We  are  now  ready  to  integrate  by  parts  in  (27 )♦  There  is  no  doubt 
that  we  can  do  this,  since  &o(t)  is  bounded  and  since  the  integrand  can  be  written 
as  the  sum  of  a  finite  number  of  tenns  in  the  form 


(32)      (1  +  t)-°  -e-^* 


k=0  \   '^   y 


0^t)  . 


Observe,   incidentally,   that  for  n  >  1  the  integrand  (32)  vanishes  at  t  ■  0  and 
t  ■   00.     For  n  >  1  we  have,   therefore. 


00 


R     »  c  / 


o^itls^^    (c.lc-l^j^k^,),. 


(c  +  n  -  l\     /  00 


ii^^  0"-^(t)dt. 


Therefore, 

a 
(33)   |R„|  <  |c|My 


-at      00 


1+t 


k"n 


/c  +  k  -  IJ 


0'^(t)dt  +  M 


n 


(c  +  n  -  11 

n 


y^-at^i. 


t)dt. 


Since  the  terms  in  the  infinite  series  are  greater  than  or  equal  to  0,   we  may  invert 
the  order  of  summation  and  integration  in  (33)  to  obtain 


(3U)  |Rj   <M|c|S^  +  MT^, 


where 


(35) 


(36) 


00 

k=n 


{'l-)/'^^^^. 


fc  +  n  -  H 


00 


e"^V-l(t)  dt  . 


To  complete  the  proof  of  the  theorem,   it  is  sufficient  to  show  that  S    -<>  0  and 

T  -e»  0  as  n  -^  00. 
n 


00 


First  consider  S  •  Our  task  is  to  show  that  S  is  the  remainder, 
n  n  ' 

,  of  a  convergent  series.  For  k  >  2a,  we  write 
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/ 


(37) 

We  have,  for  k  >  2a, 


-at  ^k^^^  dt 


Tft 


0^(log  -r-)  =  exp^  k  log  0  (log  ~) 


2n' 


(38) 


exp 


<  exp  1 


k  log  [l  -  (1  +  lof 


-  (1  +  log  ^). 2a  I 


2a^  k-ir 


■  aAk 


-2a 


where  aA  is  a  finite  constant  depending  upon  a,  but  not  upon  k.  By  the  steadily^ 
increasing  property  of  ^(t),  the  first  integral  on  the  right-hand  side  (37)  is, 
therefore,  less  than 


00 

(39)         0"(log  -^)  /  e"^"^  dt  <  Ak"^   . 


0^1og  ^)  f 
o 


For  the  second  integral  on  the  right-hand  side  in  (37)  we  have  the  upper  bound 


00 


e   -<1  -  (1  +  log  ■^)  e 


00 


_1    / 

log  It  <..iL 


-at  , 


dt 
ITt 


1  -  (1  .  log  ^)  e-^ 


L 


k 

y  dt 


2a   log  ^ 

After  the  change  of  variable  'V-  (1  +  log  ■^)  e"",  we  see  that  the  last  expression 
is  less  than 


k  N  -t 


(UO) 


<-«^'""A"-^<--^^--R£r'-^^'" 


-l-a 


Using  the  estJmates  (39),  (UO)  for  the  integrals  in  (37),  we  find,  from  the  de- 
finition (35), 
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(Ul) 


00 

s  <  A  r" 


+  k  -  i*) 


,-2a 


00 


k"n 


ti") 


Ij£lQ^  (log  ^) 

r (a+k+l)       2^ 


-1-a 


Since,  according  to  a  well-known  property  of  the  1  -  function. 


/c  ♦  k  -  1^  .  r(c  * 

\    ^     /   r(c)r 


k)    -  k 


c-1 


r(k+i)      r(c) 


as  k  ->  00, 


there  exists  a  constant  B,  depending  upon  c  but  not  upon  k,  such  that 


(U2) 


/c  +  k  -  l\ 


<  B  k 


a-l 


Similarly,  there  is  a  constant  B  such  that 


(U3) 


r(a)  Pdcn)  ^  g'  j^-a 
r(a+k+l) 


Using   (U2)  and  {h3)  in  (Ul),   we  find 


00 


-1-a 


•  ^  .-1. 


-1-a 


S    <  AB  r"    k"-^"*"  +  BB    F  k"-^(log  i-) 
k"n  k-n 


Since  both  of  these  series  have  finite  sums,  we  have 


(hh)  S  -!>  0   as   n  ->  00  . 

n 


Now  we  have  only  to  show  that  T  -oOasn-t>oo,  We  have,  for  any 


number  K  >  0, 


T  -  n 
n 


*"")|(/y")'""'" 


(t)  dt 


00 


<  B  n'^  U"^^{K)   .  i  ♦  /  e-"Ml  -  (l+K)*"* 


n-1 
?    dt 
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Since  5f(K)  <  1,  it  follows  that 


T^    <     0(1)  +  b"      /  e~"*  U  -   (1  ♦  K)  ^'\        dt- 

Setting  T  -  (1  ♦  K)  e"*,  we  find 

1 

Tn  <  0(1)  ♦  ^n^l*'^)"^/  t"-^(l-t)°-ldt 

0 

-  0(1)  +  B^d+K)-*^   ^^^)  ^^"^ 

""    r(a+n) 

-  0(1)  +  B*(l+K)'^»  B^n-l)"*"; 

here  B  and  b'  have  the  same  meaning  as  in  (U2)  and  (U3).  Letting  n  -i>  oo, 

we  find 

lim  sup  T^  <  Bb'(1+K)~** 
n  -ooo    " 

Since  K  is  an  arbitrary  positive  number,  the  right-hand  side  of  this  inequality 
may  be  made  arbitrarily  small.  Therefore, 


T   ->  0    as   n  ->  00. 
n 


This  completes  the  proof  of  the  theorem. 
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